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We show that the Wiedemann-Franz law, which holds for Landau Fermi liquids, breaks down in 
a large-n treatment of the t-J model. The calculated ratio of the in-plane thermal and electrical 
conductivities agrees quantitatively with experiments on the normal state of the electron-doped 
Pr2-iCei!Cu04 {x — 0.15) cuprate superconductor. The violation of the Wiedemann-Franz law in 
the uniform phase contrasts with other properties of the phase that are Fermi liquid like. 
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A recent experiment that measured the electrical and 
thermal conductivities of the copper-oxide superconduc- 
tor Pr2_a;Cea;Cu04 in its normal state found striking de- 
viations from the Wiedemann-Franz law^. Simply stated, 
the Wiedemann-Franz law says that fermion quasiparti- 
cles transport both electrical and heat currents, with the 
ratio of the heat conductivity k to the electrical conduc- 
tivity a given by: 




Ordinary Landau Fermi liquids respect the Wiedemann- 
Franz law, so deviations from it indicate the presence of 
non- Fermi Hquid physics. The results of Ref. 1 are there- 
fore broadly consistent with other experimental evidence 
that points to non-Fermi liquid behavior in the cuprate 
phase diagram. 

One possible interpretation of the breakdown of the 
Wiedemann-Franz law is that the quasiparticle fraction- 
alizes into separate carriers of charge and spin. To see 
what effects such fractionalization might have, consider 
first heat transport in a system of weakly-interacting elec- 
trons, which could be viewed as a crude approximation 
to electrons in the highly overdoped region of the cuprate 
phase diagram. The electrons transport both charge and 
heat. The electrical conductivity is given by the Drude 
formula a — ne'^r/m and the thermal conductivity k (x T 
because while each quasiparticle carries fixed charge e, it 
only carries an energy of order the temperature. Next 
consider a model for the undoped cuprates: the Neel 
ordered antiferromagnetic insulator with zero electrical 
conductivity. Phonons and spinwaves transport heat, 
and as both excitations are bosonic in character with lin- 
ear dispersions at low energy, each contributes similarly, 
yielding a thermal conductivity k ex. T^. 

Leaving aside these ordinary states of matter, consider 
the case in which the spins in the insulator, instead of or- 
dering, fractionalize into spinous with an extended Fermi 
surface and a non-zero density of states^'^. Now spins 
transport heat in much the same way as charges do in 
the non-interacting electron system, with k cx T. The 
Lorenz ratio is infinite, and remains significantly larger 
than one upon doping with holes or electrons. 

Just this scenario is predicted in a large-n treatment 



of the t-J model on the square lattice. We follow the ap- 
proach of Refs. 4,5,6. Lnplementing the single-occupancy 
constraint by introducing slave-boson operators bi , the t- 
J model may be written: 

H = -t [cl'^hc.Jj + H.c.) 

<i,J> 

+ \Y. ^(l^^'-^ll) ■ (2) 

Number and spin operators are related to the electron 
creation and annihilation operators by rii = c\°'cia and 
Si = (l/2)c|"af where there is an implicit sum over 
pairs of raised and lowered Greek indices. The single- 
occupancy constraint is now holonomic, blbi + cl"cia = 1, 
with the physical meaning that only a hole, or a single 
electron, may occupy each site of the lattice. (In the 
following we generally refer to hole doping with the un- 
derstanding that our calculations apply equally well to 
electron doped systems.) Included in Eq. 2 is the off- 
site Coulomb repulsion V{r); however, in the uniform 
and staggered flux phases discussed below it plays no 
role other than to contribute an additive constant to the 
total energy. 

Because spin-exchange involves no net flow of charge, 
electrical current only arises from the hopping term. The 
continuity equation relates the time rate of change of the 
occupancy on a given site to the lattice divergence of the 
current flowing into the site: 

e*^!^ = _j*e*K, H] = -y -^^-^+' ~ -^^-'-^ (3) 
dt L J. J ^ ^ w 

e—x,y 

where a is the lattice spacing between copper atoms. 
Thus 

jlj+e = -ie*t*a* {cl°'bjCj+eab]+^ - H.c.) (4) 

is the electric current flowing from site j into site j + e 
along the link connecting the two sites. We emphasize 
that neither the spin-spin exchange interaction J nor the 
Coulomb interaction V{r) appear in the expression for 



2 



the electrical current. This result, which is a direct con- 
sequence of the gauge invariancc of the spin-exchange 
and Coulomb interactions in our microscopic calculation, 
contrasts with that obtained recently^ within the more 
phenomenological "d-density wave" picture^. 

The heat current can be found by taking the 

time derivative of the Hamiltonian density: 



dt 



E 



7"^ . 



(5) 



where 



hj= J2 i-t cl"biCj^b] + H.c.) + J Si ■ Sj (6) 

i=j±x,y 

and the sum is only over the four sites i that are nearest- 
neighbors of site j. We have dropped the Tiirij interac- 
tion terms. These do not contribute to the DC thermal 
conductivity because (t?,;) remains unchanged in the pres- 
ence of currents. However, in contrast to the electrical 
current, the heat current has contributions both from 
hopping, and from spin-exchange: 



jlj+g = t*a*[c]"'bjdticj+sc.b]^g) + H.c.] 



(7) 



The model generalizes from the physical case of n = 2 
(up and down spins) to arbitrary (even integer) values 
of n by letting the Greek indices run over 1, • • • ,n. In 
the large-n limit the spin-spin interaction factorizes in 
the particle-hole channel. Formally this factorization is 
implemented via a Hubbard-Stratonovich transformation 
within the functional integral approach. Complex- valued 
mean-fields along the bonds are introduced: 



J I '{a \ 

Xij — n '^^ 



(8) 



The Xij fields function as the order parameter, and 
as they are spin-rotation invariant, there is no possi- 
bility of Neel or other spin order, and the mean-field 
Hamiltonian respects spin-rotational symmetry. Further- 
more, at sufficiently low temperatures the (holon) bo- 
son fields condense^, and we may make the replacement 
bi = b\ = y/x where x is the hole density. The mean-field 
Hamiltonian, which is exact in the n — > 00 limit, then 
takes the form: 

Hmf= ['^\Xio?-{tx + Xij){cfcioc + H.c.)]. (9) 



<i,J> 



For parameters appropriate to the cuprates, t = 0.44 
eV (following Hybertsen et al.^°) and J = 0.13 eV (ob- 
tained by Singh et al}^), the minimum energy config- 
uration has Xij both real and constant when the dop- 
ing exceeds x > 0.12. There are no broken symmetries 



in this "uniform" phase. Upon suppressing dimeriza- 
tion with the addition of a biquadratic spin exchange 
interaction^' a staggered flux (SF) phase with counter- 
circulating orbital currents'^'^'^'^'^ occurs for x < 0.12, 
that is, in the underdoped region of the phase diagram. 
(The biquadratic interaction which simultaneously ex- 
changes four fermions disappears in the physical n = 2 
limit of up and down spins, and thus does not alter the 
physics of the t-J model"''.) More realistic models in- 
clude non-zero next- and third-nearest neighbor hopping 
amplitudes'^ but these terms do not qualitatively affect 
the phase diagram or transport behavior. 

The expressions for the two currents simplify in the 
large-n limit. As the boson operators may be replaced 
by the c-number y/x, the electrical current becomes: 

jjj+i = -ie *t*x*a* {cYcj+ea - H.c.) . (10) 

Upon further replacing the fermion bilinear operator 
cj"cja with J Xi the heat current also simplifies'^: 

Jlj+g = {t*x + x)*a* [cfdtCj+^oc + H.c.) . (11) 

As the heat current differs only by the {t*x+x) prefactor 
from that of a non-interacting tight-binding system, in 
the low-temperature limit the Lorenz ratio for in-plane 
transport is simply: 



K TT 
^ " "3 
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(12) 



Thus for any x 7^ the ratio differs from unity, indicat- 
ing a breakdown of Fermi liquid theory. Note that all 
details of the scattering mechanisms cancel out in the 
ratio. Direct calculation of the two conductivities in lin- 
ear response shows that the integrals over momentum 
have identical form. Only the frequency integrals differ; 
hence for static impurities in the weak-scattering limit 
the Lorenz ratio is given by Eq. 12. We note that while 
the order parameter x is perturbed by the application of 
external electric fields or thermal gradients, the pertur- 
bation in X does not alter the DC response'^. 

In the SF phase the x-fi(3lds are complex nunibers*'^'^, 
and the prefactor [tx+xY should be replaced by |ta;-|-x|^- 
Specifically, for x = |xl exp(i^), with phase 9, the Lorenz 
ratio generalizes to'^: 



L 



\X\ 



2tx\x\ cos(e) 



{tx) 



(13) 



where L = K/{aT) and Lq = {n'^ /3){kB/e)'^ « 2.45 x 
10~^Wr^K~^ is the Lorenz number. 

We turn now to a comparison of our predictions, 
Eqs. 12 and 13, with experiments. In Fig. 1 we plot the 
Lorenz ratio as a fimction of the doping. As expected, 
the ratio diverges as the insulating limit a; ^ is ap- 
proached because the spinous transport only heat, not 
charge. In the opposite limit of large doping x ^ 0, 
and the ratio approaches unity. Landau Fermi liquid 
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theory is recovered in the dihite hmit of widely spaced 
electrons. We emphasize that the uniform phase with 
Xij constant and real does not break any symmetries. It 
exhibits weak pseudogap behavior because, according to 
the mean- field equations, |x| increases slightly in size at 
low temperatures, which in turn increases the quasipar- 
ticle bandwidth (see Eq. 9) and decreases the density of 
states (DOS)^^. For example, at a hole concentration 
of X = 0.15, Ixl = 0.024 cV at T = 500 K rising to 
Ixl = 0.026 eV at zero temperature; consequently the 
DOS drops by 2%. This contrasts with the strong pseu- 
dogap behavior of the SF phase which has a gap along 
most of the Fermi surface and which breaks time rever- 
sal invariance. In either phase, however, the fermionic 
quasiparticles are non-interacting in the n ^ oo limit and 
hence behave as long-lived Landau quasiparticles such as 
those found in ordinary Fermi liquids. 

We note that the Wiedemann-Franz law is strongly 
violated in s-wave superconducting states because while 
Cooper pairs carry charge, the condensate has no en- 
tropy. In a d-wave superconductor, quasiparticle excita- 
tions at the nodes result in a modified Wiedemann-Pranz 
law"^^. The violation that we find occurs in the normal 
state, and is a consequence of the spin-charge separation 
inherent in the large-n solution of the t-J model, and not 
of any incipient superconducting tendencies. 
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FIG. 1: The Lorenz ratio (Eqs. 12 and 13) as a function of the 
doping for t = 0.44 eV, J = 0.13 eV. The ratio approaches 
unity in the dilute limit, x —^ 1. 

In Fig. 2 we plot the temperature dependence of the 
Lorenz ratio for three dopings at which transport exper- 
iments have been conducted: x = 0.06, 0.15, and 0.26. 
For a single crystal of the La2-a;Sra;Cu04 material with 
hole doping x ~ 0.06, the resistivity was measured upon 
suppressing the superconductivity by application of a 18 
T magnetic field along the c-axis^". The thermal con- 
ductivity was, however, measured in the superconducting 
state, so it is not possible to extract a real Lorenz ratio. 
Nevertheless it is intriguing that L/Lq ~ 5 at low tem- 
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FIG. 2: The Lorenz ratio (Eqs. 12 and 13) as a function of 
temperature at dopings x = 0.06, 0.15, and 0.26 for t = 0.44 
eV, J = 0.13 cV. In the underdoped x = 0.06 case there is 
a transition to the staggered flux phase at a temperature of 
approximately 300 K. 



peratures, based on the numbers appearing in the inset 
to Fig. 3 of Rcf. 20. This compares reasonably well with 
the theoretical value of 4.1 seen in Fig. 2. 

At optimal doping, experimentally available magnetic 
fields can only eliminate superconductivity in electron- 
doped compounds. In Ref. 1 a 13 T field was applied to 
Pr2_a;CexCu04 at a; = 0.15 to access the normal state. 
The measured ratio of L/Lq ~ 2 found at temperatures 
above 0.3 K is again in reasonable quantitative agree- 
ment with the theoretical value of 1.95. At temperatures 
below 0.18 K, however, the experimentally determined 
ratio drops rapidly below one. We have no explanation 
for the observed behavior at the lowest temperatures^^ . 

Finally, in the highly overdoped regime Proust et al?"^ 
have studied the Tl2Ba2Cu06+5 material at a hole con- 
centration of .T » 0.26. Superconductivity was sup- 
pressed in a 13 T field, and L/Lq = 0.99 ± 0.01 in good 
agreement with the Wiedemann-Franz law for Fermi liq- 
uids. The theoretical value of the ratio is 1.5. We spec- 
ulate that the persistence of non-Fermi liquid behavior 
at large doping in the mean-field theory is an artifact 
of the large-n approximation. Finite-n corrections could 
possibly restore Fermi liquid behavior in the overdoped 
region. At large doping |x| is small compared to the ef- 
fective hole hopping amplitude tx, so fluctuations in x 
may be expected to be relatively more important than at 
low doping. 

In summary we have shown that the Wiedemann-Franz 
law is violated in a mean-field treatment of the t-J model. 

Our analysis, which holds for weak scattering, is exact in 
the n ^ oo limit. The Lorenz ratio is significantly larger 
than one both in the uniform phase {x > 0.12) and in 
the SF phase {x < 0.12). The theoretical prediction is 
in reasonably good quantitative agreement with existing 



4 



experimental measurements on the cuprate materials. 

Note added: After this work was completed a paper 
by Kim and Carbotte (KG) appeared^^ that examined 
the Wiedemann- Franz law within the context of the phe- 
nomenological d-density wave picture. There are several 
differences between their work and ours. The main dif- 
ference is that we study both the uniform phase which 
has no broken symmetries, and the SF or d-density phase 
with time-reversal breaking coimter-circulating currents. 
We find that the Lorenz ratio is significantly larger than 
one in both phases. Furthermore, at low temperatures 
KG find only small deviations from the Wiedemann- 



Franz law. This is due in part to the fact that their d- 
density order parameter (the analog of our Xij ) was cho- 
sen to be purely imaginary (equivalent to setting = 7r/2 
in our Eq. 13) and also because their kinetic energy is not 
rescaled by the slave-boson doping factor, x, as it is in 
our microscopic analysis of the t-J model. KG also find 
a large temperature variation in the Lorenz ratio for the 
case of strong scattering because the quasiparticle life- 
time has a strong frequency dependence. 
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